Riemann— Hilbert analysis for Laguerre polynomials 
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Abstract. We study the asymptotic behavior of Laguerre polynomials L^ n \nz) as n — » oo, where 
a n is a sequence of negative parameters such that —a n /n tends to a limit A > 1 as n — > oo. These 
fS) 1 polynomials satisfy a non-hermitian orthogonality on certain contours in the complex plane. This 
fact allows the formulation of a Riemann-Hilbert problem whose solution is given in terms of these 
Laguerre polynomials. The asymptotic analysis of the Riemann-Hilbert problem is carried out by 
the steepest descent method of Deift and Zhou, in the same spirit as done by Deift et al. for the 
Qh' case of orthogonal polynomials on the real line. A main feature of the present paper is the choice 
[ of the correct contour. 
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*S ■ 1 Introduction 

1.1 Generalized Laguerre polynomials 

The classical Laguerre polynomials Ln are orthogonal on the interval [0, oo) with respect 
OO ' to the weight x a e~ x , that is 
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poo 

/ L£Xx)L%\x)x a e- x dt = Q, lin^m. (1.1) 
Jo 



(N ' 

^t- : 
o . 

The integral in (1.1) converges only if a > — 1. The Laguerre polynomials are given by the 



explicit formula 



fc=0 

and by the Rodrigues formula 



u—n \ / 



L^\x) = ±aTV (f] n [x<* +n e~*] , (1.3) 



n\ \dx y 

see e.g. [P0| . Both ( |1.2| ) and ( |1.3| ) make sense for arbitrary a (even complex) and they define 
the generalized non-classical Laguerre polynomials. The recurrence relation 

-xLl a \x) = (n + l)L { :l(x) -(2n + a + 1)4%) + (n + c^L^x), (1.4) 

with Lq = 1, L_l = 0, and the second order differential equation 

xy"{x) + (a + 1 -x)y'(x) +ny(x) = 0, y(x) = L^(x), (1.5) 

continue to hold for arbitrary a. We consider in this paper only real and negative a, although 
extensions to complex a are possible. 
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1 



1.2 Relations to other polynomials 



Laguerre polynomials with negative parameters appear in the literature in a number of forms. 
First we note the special cases a = —n and a = —n — 1, where we have 



nl 



;i.6) 



and 



k=0 



(1.7) 



respectively. Thus for a = —n — 1, the generalized Laguerre polynomial agrees (up to a sign 
if n is odd) with the partial sum of the exponential series. More generally, we have that for 
a = —n — m with n,m G N, the generalized Laguerre polynomials appear as the numerator 
and denominator polynomials in the rational Pade approximant for the exponential function. 
To be precise, if 



Pn,m{x) = (-1)' 



and 



n + m 
n 



n + m 



m 



then p n>m (0) = q n ,m(0) = 1, and 

Pn, m (x) - q n , m (x)e x = O (x n+m+1 ) as x -> 0, 



;i.9) 



;i.io) 



see e.g. p3| , 

Laguerre polynomials with negative parameters are also related to the so-called general- 
ized Bessel polynomials 



y n (z; a) = 53(?) 



[n + a — 1) 



since 



y n (z;a) = (-lTn\ (|)" L^~ a ^ 



1.121 



The usual Bessel polynomials correspond to a = 2 in ( l.llj) . See Grosswald |T^| for a 
comprehensive account of these polynomials. 



1.3 Earlier work on asymptotics 

For a > — 1, the Laguerre polynomials satisfy the orthogonality (|1 . 1| ) on the positive real 
axis, and therefore they have only positive real zeros. This property is lost for a < — 1. 
Indeed, the generalized Laguerre polynomials may have many non-real zeros. In Figure [I] 

for a number of values of A > 0. Similar plots are 



we have plotted the zeros of4o 40A) (4(k) 



shown in the paper [ 20| . 
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Figure 1: Zeros of generalized Laguerre polynomials L 
(left), A = 1.01 (middle), and A = 2 (right). 



(-An) , \ 



for n = 40 and A = 0.81 



In Figure [I] we see that the zeros cluster along certain curves in the complex plane. 
Martinez et al. |20| identified these curves as trajectories of a quadratic differential, depending 
on a parameter A. For A > 1, the curve is a simple arc, which as A decreases to 1 closes 
itself to form for A = 1 the well-known Szego curve p9|, [24]. For < A < 1, the curve 



consists of a closed loop together with an interval on the positive real axis. 

A number of rigorous results on the asymptotic behavior of the zeros of generalized 
Laguerre polynomials Li? n \nx) such that 



lim — = —A 

n-^oo fl 



;i.i3) 



are known from the literature. The first result is due to Szego |29[ who studied the partial 
sum of the exponential series, that is a n = —n—1 see ( |1.7|) . Szego showed that the normalized 
zeros tend to the curve which now bears his name. Olver |21| considered the zeros of Hankel 
functions, which includes the Bessel polynomials as a special case. In terms of the generalized 
Laguerre polynomials, this is the case a n = —2n — 1. Saff and Varga P6| studied the zeros 
and poles of Pade approximants to the exponential function, see (|1.8| )- (|1.9| ). Their main 
result says that for integers a n < —n such that ( [1.1 3D holds, all zeros of L { n n \nx) tend to 
a well-defined curve as n — > oo. The curve depends on A > 1 only, and coincides with the 
curve described in ||20|| . Saff and Varga also obtained the weak limit of the zero counting 
measures. The proofs in [^6| can be extended without any difficulty to non-integer a n < —n. 

Stated in terms of generalized Bessel polynomials 
zeros are due to De Bruin et al. ||, Carpenter 
papers deal with the limit ( |1.13|) with A = 2. The latter paper also presents uniform 
asymptotic expansions of the generalized Bessel polynomials. In very recent work, Dunster 



1.11 



1.12|) asymptotic results on 

These 



and Wong and Zhang BT 



II | establishes uniform asymptotic expansions in the complex plane for the case of general 

and A 



A, with the exception of A 
case A > 1. 



1. The results on zeros in 14 are restricted to the 



1.4 Asymptotics from Riemann— Hilbert problems 



Most papers cited above use some form of the steepest descent technique for integrals, see 
especially [^6| and Martinez et al. |20[ use an orthogonality relation in the complex 



plane satisfied by generalized Laguerre polynomials. The approach of Dunster |14| starts 
from the differential equation ( Oj ) and is based on techniques developed by Olver |22 . 

In this paper we derive asymptotics of generalized Laguerre polynomials using the non- 
linear steepest descent / stationary phase method for Riemann-Hilbert problems introduced 
by Deift and Zhou in [0, and further developed in |13] and [11|. In later developments, 



the method was applied successfully to problems in random matrix theory ||, ||, and in 
orthogonal polynomials ||, ||, ||10|| , fT9|| , and combinatorics ||. For review of some of 
these developments, and a pedagogic introduction to some of the material of random matrix 
theory, orthogonal polynomials, and Riemann-Hilbert problems, see [BJ. 

The Riemann-Hilbert approach to the asymptotics of generalized Laguerre polynomials 
starts from the observation that these polynomials satisfy orthogonality relations in the 
complex plane. The orthogonality is on a contour E going around the positive real axis, 
but otherwise being quite arbitrary. The orthogonality property allows the formulation of 



a Riemann-Hilbert problem, due to Fokas, Its, and Kitaev whose solution is given 



in terms of Ln '■ The Riemann-Hilbert problem is analyzed in the large n limit with the 
steepest descent method as done in M and [TIJ for orthogonal polynomials on the real line. 

A novel feature for the problem at hand is that the arbitrary contour E has to be chosen 
in a correct way in order to arrive at a Riemann-Hilbert problem which is amenable to 



subsequent asymptotic analysis. The correct contour was described in [20]. It is a curve 



with the S-property of Stahl and Gonchar and Rakhmanov [I5|. The structure of the 
curve depends on the value of A as already explained before. In this paper we analyze the 
case of an open contour, that is, the case A > 1. In subsequent work we consider the case of 
a closed loop plus an interval (0 < A < 1) and the case of a single closed contour (A = 1). 
We note that the steepest descent / stationary phase method for Riemann-Hilbert problems 
was augmented to handle cases in which the contour selection involves determining a set of 



nontrivial curves in the plane in fL8|| , by Kamvissis, McLaughlin, and Miller, in the context 
of the semi-classical limit of the focusing nonlinear Schrodinger equation. 

We emphasize that the main interest in the present paper lies in the method we use and 
not in the results obtained for the Laguerre polynomials. In particular, we do not improve 
upon the asymptotic expansions of Dunster |14[]. The steepest descent method for Riemann- 
Hilbert problems is a very powerful new method, and its use in the study of classical special 
functions is new. In future work we consider generalized Laguerre polynomials for the cases 
< A < 1 and A = 1, and the Riemann-Hilbert approach will lead to new results for these 
cases. 



2 Complex orthogonality and the formulation of the 
Riemann-Hilbert problem 

2.1 Orthogonality 

For a < — 1, the generalized Laguerre polynomial Ln is not orthogonal on the positive real 
axis, but instead satisfies a non-hermitian orthogonality in the complex plane. 

Let T be the collection of all simple Jordan curves E in C \ [0, oo) that are symmetric 
with respect to the real axis, and such that there is M > such that for all x > M, there is 
y(x) > 0, such that the intersection of E with Kez = x consists of the two points x ± iy(x), 
and lim^oo y(x) = L exists and is finite (possibly 0). Any curve E e JF divides the complex 
plane into two domains, Q + and where f2_ contains the positive real axis. We choose 
the orientation of E such that f2+ is on the +- side (i.e., on the left) while traversing E and 

is on the — side. So E is oriented clockwise as in Figure ||. 

In what follows we define x a with a branch cut along the positive real axis. Thus x a = 
|x| a e margx ' with argx E [0,2%). 
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Figure 2: Example of contour £ 
Lemma 2.1 LetY, <E n <E N, and «6l TTien 

y 4 a) (aOxVe - * = 0, /or A: = 0, 1, ... , 
If in addition a + n + 1 N, i/ien 

y L^^x^e^dx^O, fork = n. 



n — 1. 



(2.1) 



(2.2) 



Proof. The orthogonality fl2.1| ) follows from the Rodrigues formula ( |1.3| ) by repeated in- 
tegration by parts. In the same way, we also get 

y L^\x)x n x a e- X dx = ^ y (x a+n e- x )x n dx={-l) n J x a+n e~ x dx 

For a + n > — 1, we deform £ to the positive real axis to obtain 



p poo 

/ L^\x)x n x a e- X dx = (-l) n (l-e 2 ™) / x a+n e~ x dx 
Jt, Jo 



l) n+1 2ie ma sin(7ra)r(a + n + 1), 



(2.3) 



where V denotes the Gamma function. By analytic continuation the integral in fl2.2|) is equal 
to ( |2.3| ) for every a, and ( |2.2| ) follows. □ 



The formula (|2.1|) expresses orthogonality with respect to the complex measure x a e x dx on 
E. 



2.2 Riemann— Hilbert problem 

Let a G R. We consider the monic polynomials 

'-l) n n\ 



Pn{z) 



n" 



n = 0,1, 



(2.4) 



Introducing a change of variables x = nz in fl2.1p and Q2.2| ), we see that 

f = 0, for Jfc = 0, 1,... ,n- 1, 

P„(z)z fc z a e- n2 ^<^ (2.5) 
s I 7^ 0, for k — n, 

for every contour E G JF, provided that a + n + 1 ^ N. 

The polynomial P n is characterized through a Riemann-Hilbert problem due to Fokas, 
Its, and Kitaev |T5[ . 



Riemann-Hilbert problem for Y: 

Let E be a contour from the class J 7 , that divides the complex plane into two parts Q + and 
fi_, as above. The problem is to determine a 2 x 2 matrix valued function Y : C \ E — > C 2x2 
such that the following hold. 

(a) K(^) is analytic for z G C \ E, 

(b) K(z) possesses continuous boundary values for z G E, denoted by Y + (z) and K-(^), 
where and Y_(z) denote the limiting values of Y(z') as z' approaches z G E from 
Q + and fi_, respectively, and 

/ I z a e~ nz \ 

Y + (z) = Y.(z)l L Q \ J, forzGE, (2.6) 

(c) Y(z) has the following behavior as z — >• oo: 
= fj + O fi^ f asz^oo, zG<C\E. (2.7) 



Proposition 2.2 Lei a G K with a + n G" N. JTien i/ie unique solution of the Riemann- 
Hilbert problem for Y is given by 



Y(z) 



( Pn(z) 

V 



- f 

2m J 



a -nC, 



d( 



(2- 



where P n (z) is the monic generalized Laguerre polynomial 

\n+l rn n+a 



/ 

and 



Qn-i(z) 



(-l) n+1 Jl n+a TC-™ r(a) 



sin(7ra)r(a + n) 



L y n l x {nz) 



(2.9) 



Proof. The proof is as in Section 3.2]. Here we will only give the proof for the second 
row, since that is where the condition on a plays a role. 

From (|2.6|) it follows that Y 2 \ is an entire function, which by ( |2.7|) satisfies Y 2 \{z) = 



0(z 



as z 



oo. Therefore Y 2 \{z) = Q n -i(z) for some polynomial Q n -\ of degree at most 



n-l. The (2,2) entry of the condition (£§) is (F 22 )+(^) = (Y 22 )-(z) + Q n - l (z)z a e- nz , which 
by the Sokhotskii-Plemelj formula yields 



Yoo(z) 



1 

2vri 



Qn-i(C)C Q e 



a -T»C 



d(. 



(2.10) 



From the (2,2) entry of Q2.7p it follows that Y 2 2(z) = z n + 0{z n x ) as z — > oo. Because of 
( p.lO| ) this gives the conditions 



and 



J Q n ^(C)Ce- n( C k d( = for A; = 0, . . . , n - 2, 

^g n _ 1 (C)C Q e-<c n - 1 ^C = -2«. 



(2.11) 



(2.12) 



The orthogonality conditions (|2.11|) are satisfied if Q n -\(z) = c? n L„_ 1 (n2;), and the constant 
d n must be chosen so that ( |2.12| ) holds as well. Because of Lemma 2.1 this can be done if 
a + n £ N, and the result is given by (|2.9|) . □ 



3 Selection of the right contour 
3.1 The right contour E 

In Section 3-6, we assume n G N and a < —n are fixed. We write A = —a/n, so that A > 1. 
[Later, when we let n —>■ oo, a and A, as well as other notions introduced in these sections, 
will depend on n] 

A major step in the analysis of the Riemann-Hilbert problem for Y is the selection of the 
right contour. In order that the subsequent analysis works, the contour cannot be arbitrary 
but has to chosen in a precise way. The contour depends on A. We define 



f3 = 2 - A + 2iVA - 1 



(3.1) 



Martinez et al. 201 showed that the values of z for which 



1 

2vri 



0) 



1/2/ 



0) 



1/2 



ds is real, 



where the branch of the square roots is chosen so that they are analytic and single valued 
on the path of integration from (3 to z, form a system of curves as shown in Figure [5] for a 
number of values of A. In geometric function theory these curves are known as trajectories 
of the quadratic differential iS^m^zH (p s ^ see |pg|| . We see four smooth (in fact analytic) 





Figure 3: Curves where ^ f| ^f^- ds is real, for the values A = 1.01 (left) and A = 2 (right). 
The two points of intersection are j3 and j3. 

curves. Two curves are connecting (3 with (3, one of them crosses the negative real axis, and 
the other one crosses the positive real axis. 
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Definition 3.1 We define T as the trajectory of the quadratic differential - $ d 2 s from 
P to (3 which crosses the negative real axis. V is oriented from (3 to (3. 

We put 

R(z) = {z-f3) 1/2 {z-P) l/ \ zeC\T, (3.2) 

where the branch is chosen which is defined and analytic on C \ V, and which is such that 
R(z) ~ z as z — > oo. For s G T, we use R+(s) and R-(s) to denote the limits from the 
+-sides and — sides, respectively. As usual, the +-side of an oriented curve lies to the left, 
and the — side lies to the right, if one traverses the curve. 
Then by definition of T, we have 

J_/**±Md. is real for every , eV, (3.3) 
2m Jp s 

with integration along the +-side of V . 

It is also of interest to know where 7^7 Jjf ds or ^7 fp ds is purely imaginary. 
These are the dotted curves shown in Figure [|. The dotted curves are analytic extensions 
of the solid ones. 



Figure 4: Curves where 7^7 ds is real (solid lines) and curves where ^7 Jb ds or 

J/3 ~r~ ds is purely imaginary (dotted lines), for the values A = 1.01 (left) and A = 2 
(right). 

We can now state which contour S to choose. 

Definition 3.2 We let E be the contour in T consisting of T together with the two dotted 
curves that form the analytic extension of T. 

We denote the part of S \ T in the lower half-plane by and its mirror image in the 
upper half-plane by S 2 . 

So we have a disjoint union S = r U Si U S 2 . Figure |5| shows the curve E for two values of 
A, together with the zeros of the corresponding Laguerre polynomial L\ i0A \40x) of degree 
40. The figure shows that the zeros are close to T, and that they are in the domain 
These findings will be confirmed by our final result, Corollary 7.2 below. 

Remark 3.3 We have chosen E 2 so that J* ^j^- ds is real and positive on E 2 . This is not 
essential. What is important for the subsequent analysis is that it has positive real part on 
E 2 . This means that we have the freedom to deform E 2 , as long as we take care that the 
real part of ds is positve on E 2 . However, it will be convenient to do this deformation 

only away from (3, so that in a neighborhood of (3, we have E 2 exactly as we defined it. 
Similar remarks apply to ^j^- ds and Ei. 



x 



Figure 5: The curves S = F U £1 U E 2 , and the zeros of L^q 40 ^ (40x), for the values A = 1.01 
(left) and A = 2 (right). 

3.2 A probability measure on T 

The following proposition gives one of the crucial properties of T. 

Proposition 3.4 The (complex) measure _L R +( s ) d s ^ s a probability measure on V. 
Proof. We show first that 

1 rRM d s=l. (3.4) 



2iri J r s 

Let I = J r R+ j s ^ ds. Since -R-(s) = — R + (s) for s G T, we have 

27= / M-^£>U=<^ rfs . 



where 7 is a closed contour encircling the curve T once in the clockwise direction and not 
encircling z = 0. 

After contour deformation, we pick up residues at z = and at z = 00, namely 

2/ = 27uRes^-27rzRes^. (3.5) 

z=0 Z z=oo z 

The residue at z = is 

Res ^ = R(0) = \p\ = y/(2 - A) 2 + A(A - 1) = A (3.6) 

2=0 z 

and the residue at z = 00 is the coefficient of z~ x in the Laurent expansion of R(z)/z 
^ = (l-/3/,)^(l-^)V2 



Thus 



!--±w). 



Res ^(fl = -P±l = -Ref3 = -(2 - A). (3.7) 

2=00 2 



Hence by (^-Q we have 21 = 2ni(A + (2 — A)) = Am, so that Q) follows. 



Having ( |3.4| ) we can now prove the proposition. Let t i— > z(t) for i G [0, t ], be the arc 
length parametrization of T starting at (3. Thus #(0) = (3 and z(to) — P- Then we have for 
z = z(t) with £ G (0,t ), 

1 rw'i,, ' rm 1r)Jr , (3. 8) 



27TZ S 27T2 7q z[t 

and this is real for every t by construction of T. It has the value for t = and the value 1 
for t — to (due to (|3.4|)). The derivative of (|3.8|) with respect to t is i ^+^W) ^^ which is 



2ni z(t) 

not zero for t G (0, to). Thus ( [3.8|) can only increase from to 1 as t increases from to to- 
Hence ds is a positive measure on T. It is a probability measure because of (p.4|). □ 

3.3 Auxiliary functions 

With the measure d/i(s) = ^^j^-ds on T we define the so-called ^-function as follows. 
Definition 3.5 The ^-function is the complex logarithmic potential of \x, that is, 

g{z) = Jlog{z-s)dfi{s), *eC\(ru£i), ( 3 - 9 ) 

where for each s we view log (z — s) as an analytic function of the variable z, with branch 
cut emanating from z — s. The cut is taken along V U Si. 

We need two more functions. 

Definition 3.6 The 0-function is defined as 

, , in r{ s 



f ^^ds, 2GC\(ruS 1 U[0,oo)), (3.10) 

J/3 s 

where the path of integration from (3 to z lies entirely in C \ (T U S x U [0, oo)), except for 
the initial point /3. 

The 0-function is defined as 

fcz) = I I — ds, 2GC\(ruE 2 U[0, oo)), (3.11) 

2 J /3 s 

where the path of integration from /3 to z lies entirely in C \ (T U S2 U [0, 00)), except for 
the initial point (5. 

It is immediate from ( |3.3| ) that <fi+(z) is purely imaginary for z G V . By Proposition 3.4 
its imaginary part increases from to 7r as z traverses the curve V from {3 to (3. In particular 
we have (f>+(/3) = m. Similarly 4>-(/3) = —ni. Thus we have 

{(f)(z) + ni for z G 
(3.12) 
(j)(z) — ni for z G f2_. 

Proposition 3.7 There is a constant £ such that 

g{z) = ^(Alogz + z + £)~(f ) {z), z G C \ (r U E x U [0, 00)), (3.13) 
where logz is defined with a branch cut along [0, 00). 



in 



Proof. We note that 



1 

2i\i J T z — s s 



1 RJs) 1 1 
as 



1 R(s) 



2 2ixi L z — s s 



ds 



where 7 is a closed contour in C \ T, that encircles T once in the clockwise direction, but 
does not encircle z and 0. Then exactly as in the proof of Proposition 3.4, we deform the 
contour, and now pick up residues at z, and 00. We find for z G C \ I, 



9'M 



Res (J-M) + Res (J-M) _ Res (J-M 

s=z \Z — S S J s=0 \Z — S S J s=oo \Z — S S 



+ - + 1 



2; z 



(3.14) 



Let z G C \ (r U Si U [0, 00)). Integrating ( [3.14j ) from /3 to z along a curve in C \ (r U 
Si U [0, 00), we find 



g{z) = g(P) + ^(Alogz + z) - \{A\og(3 + (3) - <f>{z). 



(3.15) 
□ 



Thus Q3TT3|) holds with t = 2g{(5) - (A log/? + (3). 
3.4 Jump properties of g 

From Proposition 3.7 we obtain the following jump relations for g across the contour S. 
These jumps are crucial for the subsequent analysis. 

Proposition 3.8 (a) We have 

9+(z) - 9-( z ) = 2ni f° r z G S i' ( 3 - 16 ) 

and 

g+(z) - g-(z) = -cj> + {z) + <f>4z) = -2</» + (z) = 20_(z) for z G T. (3.17) 



(b) We have, with the same constant t as in Proposition 3.7, 

g+(z) + g~(z) = Alogz + z + £ forzeT, 



(3.18) 



and 



g + (z)+g_(z) = Alogz + z + £-2(f>(z) for z G S x , 



g + (z) + g-(z) = Alog z + z + £ - 2cf)(z) for z G S 2 . 



(3.19) 



(3.20) 



Proof. In (|3.13| ) we let z approach S, from the +- and — sides, respectively, to obtain 

g±(z) = l(Alogz + z) + X -i - 4>±(z), for z G S. (3.21) 



1 1 



Since changes sign across T, Q3.17 ) and fl3.18| ) immediately follow from (|3.21 ). For z G Si, 
we have by ( ETT2D and (|OTD 



= -0+O)+0-(>) = -(0+(z)-7ri) + (0_(z)+7ri) 

= 2tTZ-0 + (z)+0_(z), 



which is ( 3.16Q as is analytic across Si. In addition, we have 



9+{z)+g-{z) = Alogz + z + £- <f) + {z) -(j)-(z) 

= Alogz + z + i - 4>+(z) - 4>~{z) for z G S x . 
which yields flSTISD . Similarly, (|3~20D follows. □ 

4 First two transformations: F^[/^T 
4.1 First transformation Y i— ► [/ 

With the g-function and the constant £ from Proposition 3.7, we perform the first transfor- 
mation of the Riemann-Hilbert problem. 

Definition 4.1 We define for z G C \ S, 

U(z) = e ~n(m)-3 Y (z)e- n9{z)a:i e nW2)r7 \ (4.1) 

Here, and in what follows, 03 denotes the Pauli matrix °3 = ( q ^ I > so that f° r example 



e -ng(z)a 3 



e ~ng(z) Q 

e ™f0) 



From the Riemann-Hilbert problem for K it follows by a straightforward calculation that 
U is the unique solution of the following Riemann-Hilbert problem. 

Riemann-Hilbert problem for U : 

The problem is to determine a 2 x 2 matrix valued function [/ : C \ S — >C 2x2 such that 

(a) [/ (z) is analytic for 2 G C \ S, 

(b) t/(z) possesses continuous boundary values for z G S, denoted by and U-(z), 
and 

/ p -"(S+(«)-9-(2)) 7 ~An -nz n(g + (z)+g-(z)-£) \ 

for 2 G S, 

(c) £/ (z) behaves like the identity at infinity: 

f/(z)=/ + oQ as z ^00, zeC\E. (4.3) 
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The jump relation ( |4.2[ ) for U has a different form on the three parts T, Si and £2. On 
T we have that g + — g_ = — 2<p + = 20_ by (|3.17|) , and g + + g_ = Alogz + z + £ by Q3.18Q 
so that 



On Si we use ( [3.161 ) and Q3.19Q to obtain 

U+(z) = U-(z 

and similarly it follows that 

U + (z)=U(z 



,2n4> + (z) 







3 2n0_ (z) 



for z e r. 



Y e -2n<j>(z) 
1 



for z G Ei, 



1 g-2nfli(z) 

1 



for zeE 5 



(4.4) 



(4.5) 



(4.6) 



The transformation Y \— ► C/ has the effect of normalizing the Riemann-Hilbert problem 
at infinity. In addition, by the construction of E, we have that is real and positive on E 2 , 
and (f) is real and positive on Ei. So the jump matrices for U in (|4.5| ) and (|4.6| ) are close to 
the identity if n is large. Since (f) has purely imaginary boundary values on both sides of T, 
the jump matrix for U on T in (|4.4j) has oscillatory diagonal entries. 



4.2 Second transformation U 1— > T 



The jump matrix for [/ on T, see 

; 2n0+(z) J 







3 2n0_ (z) 



factors as 



1 

3 2n</>_(z) 



1 
-1 



1 

e 2n4> + {z) j 



(4.7) 



Observe that the first matrix in the right-hand side of ( |4.7|) can be analytically continued to 
the — side of the contour T, and in doing so, the (1,2) entry becomes exponentially decaying 
in n. Similarly, the third matrix in the right-hand side of (|4.7|) can be analytically continued 
to the +- side of the contour T, and in doing so, the (1,2) entry also becomes exponentially 
decaying in n. We are thus led to introduce the following "contour augmentation" step, 
as part of the steepest descent / stationary phase method for Riemann-Hilbert problems 
developed by Deift and Zhou. The oriented contour E T consists of E plus two simple curves 
£3 and £4 from (5 to (3, contained in Q + and respectively, as shown in Figure |6|. We 
choose £3 and £ 4 such that Re <f>(z) < on £ 3 and £4. 

It is possible to choose such curves. Indeed, is positive on £ 2 , and its real part vanishes 
on r and on the other solid lines shown in Figure [|. So Re > in the full region on the 
right. In the two other regions, bounded by the solid lines, we then have that Re0 < 0. 
Note that Re0 does not change sign across T. 

Note that by ( |3.12| ) we also have Re0(^) < on £ 3 and £4. 

Then C\E T has four connected components, denoted by Qi, fl 2 , ^3, and ^4 as indicated 
in Figure 0. 



Definition 4.2 We define T : C \ £ T -> C 2x2 by 

T(z) = U(z) for zGfliU 4 , 



(4.8) 
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Figure 6: Contour S T = V U Ej and the domains flj, j = 1, . . . ,4, for the Riemann- 
Hilbert problem for T, for the value A = 1.01. 

T(z) = U(z) ( _J n<j>{z) M for 2 G fia, (4.9) 

T(z) = U(z)(j <Kz) °A forzGfi 3 , (4.10) 

Then from the Riemann-Hilbert problem for U and the factorization (|4.7| ) we obtain that 
T is the unique solution of the following Riemann-Hilbert problem. 



Riemann-Hilbert problem for T: 

The problem is to determine a 2 x 2 matrix valued function T : C \ S T — > C 2x2 such that 
the following hold: 

(a) T(z) is analytic for z G C \ S T , 

(b) T(z) possesses continuous boundary values for z G S T , denoted by T + (z) and T_(z), 
and 

T+(z) = T_(z) ( _° x J) for z G T, (4.11) 
1 

e 2 ' 



T + (z) = T„(z)( ^ {z) 1 ) for^GS 3 US 4 , (4.12) 



and 



T + (z)=T4z)[ J x ) forzGEx, (4.13) 



1 p -^n<p(z) 

r+(z)=T_(z)( Q 1 ) for^GS 2 , (4.14) 



(c) T(z) behaves like the identity at infinity: 



T(z) = I + ()[-) asz^oo, zeC\E T . (J..U)) 



1 A 



5 Construction of the parametrix for T 
5.1 Parametrix away from (3 and (3 



As remarked following (|4.7| ), the jump matrices appearing in ( 4.12 ) for z G S 3 U S 4 are 
exponentially close to the identity matrix away from (3 and (3. Similarly, the jump matrices 
appearing in (|4.13|) and (|4.14j ) are also exponentially close to the identity matrix away from 



f3 and (3. This hints that these portions of the contour on which the Riemann-Hilbert 
problem for T is posed should be somehow negligible. Thus we expect that the leading 
order asymptotics is determined by the solution N of the following model Riemann-Hilbert 
problem: 

Riemann-Hilbert problem for iV: 

The problem is to determine iV : C \ Y — > C 2x2 such that the following hold. 

(a) N(z) is analytic for z G C \ T, 

(b) N(z) possesses continuous boundary values for z G T \ {/3,/3}, denoted by N + (z) and 
N_(z), and 

N+(z) = N-(z)( _° x for z G T \ {(3 J}, (5.1) 

(c) N(z) = I + O (±) for z -> oo. 

This Riemann-Hilbert problem for is solved explicitly by, see |1(| p. 1520] or || p. 200], 

a(z)+a(z) _1 a(z)—a(z)~ 

= I a ( z )-a( z )-i a(z)+a(z)- 1 ) J (5-2) 
\ -2i 2 

where 



o(2) = trfw- (5-3) 

The branches of the roots in ( |5.3|) are chosen such that a(z) is analytic on C \ T and 
lim a(z) = 1. 



Remark 5.1 The solution ( |5.2| ) is not the only solution to the Riemann-Hilbert problem 
for N. It is the unique solution that satisfies, in addition to (a), (b), and (c), the condition 

(d) Near the endpoints (3 and f3, we have 

N{z) =0(|z-/3|- 1/4 ) asz^/3, 

N{z) = 0(\z - asz^ft 
with the O-term being taken entry-wise. 



1F> 



Remark 5.2 For explicit calculations, it is useful to have an alternative expression for N. 
For z = x real, it is clear from ( |5.3|) that a(x) has modulus one. If arg(a(x)) = 6(x), then 
(|5.2|) shows 

cos6*(:r) sm9(x) 
— sinQ(x) cos9(x) 



N(x) 



Since tan(26*(x)) = — ^-2+ a > we then fi n d 



cos(| arctan( 2 ^ L -j)) — sin(| arctan(^|— |)) 
= I 2 ~ 2 _ J (5.4) 

sin (| arctan( 2 _^ +y | )) cos(| arctan( 2 _ 2 ^ )) 

first for z = x real, but then also for arbitrary z G C\T by analytic continuation. We have to 
take the appropriate branch of the multivalued arctan function in ( |5.4| ). Using trigonometric 
identities, one may then check that ( |5.4| ) reduces to 



N(z) 



V 



2 / v 2 



for z G C \ T, (5.5) 



where as usual we have R(z) — (z — /3) 1 / 2 (z — /5) 1//2 . We can also verify directly that ( |5.5| ) 
solves the Riemann-Hilbert problem for N. 

5.2 Parametrix near (3 

The next step is a local analysis around the points j3 and /3. We need to construct a local 
parametrix P in a neighborhood £7$ = {z G C | \z — j3\ < 5} of (3 such that 

• P satisfies the jumps for T exactly in Us, 

• P matches iV on the boundary of Us up to order 1/n. 

See Figure [7| for the contours S T n Ug. 
More precisely, we have 

Riemann-Hilbert problem for P: 

The problem is to determine, for a given 5 > sufficiently small, a 2 x 2 matrix valued 
function P : Z7 5 \ S r -> C 2x2 such that 

(a) -P(z) is analytic for z E Us\ S T , and continuous on Us \ S T , 

(b) P(z) possesses continuous boundary values for z G S T fl Z7<5, denoted by P+(^) and 
•P-Oz), and 

P+(z) = P-(z) ( _° x J ^ forz6rn[/ 5) (5.6) 
P+(z) = ( J ) for 2 G (S 3 U £4) fl E/j, (5.7) 



e 



and 



I p -2n<p(z) 

P+(z)=P4z)l 1 1 for^GS 2 nf/ 5 , (5.8) 



0.35 



0.25 




0.05 



0.75 0.8 0.85 0.9 0.95 



1.05 1.1 1.15 1.2 



Figure 7: Neighborhood Us of (3 and the parts of the contours T, S 2 , S3, and S 4 that are 
within Ug. The value of A is 1.01. 



(c) There exists a constant C > such that for every z G dUs \ S T , 

C 

\\P(z)N-\z) - I\\ < -, 



(5.9) 



where 



is any matrix norm. 



The construction of P follows along the same lines as given by Deift et al. [T(| . From its 
definition ( |3.1U| ) it is easy to see that the 0-function has a convergent expansion 



» = (z-/3) 3 / 2 ]Tc fc (z-/3) fc , c ^0, 



(5.10) 



fc=0 



in a neighborhood of (3. The factor [z — (3) 3 / 2 is defined with a cut along r U Si. Then 



/(*) 



2/3 



(5.11) 



is defined and analytic in a neighborhood of (3. We choose the 2/3-root with a cut along V 
and such that f(z) > for z G S 2 . Recall <fi > on S 2 . 

Then /(/?) = and /'(/3) 7^ 0. Therefore we can and do choose 5 so small that ( = f(z) 
is a one-to-one mapping from Us onto a convex neighborhood f{U$) of ( = 0. Under the 
mapping ( = f(z), we then have that S 2 fl Ug corresponds to (0, 00) n f(Ug) and that TnUg 
corresponds to (—00, 0] fl f(Ug). 

Now we specify how to choose S T near f3. For an arbitrary, but fixed a G (tt/3, 7r), we 
choose S 3 and S 4 such that / maps S 3 fl Us and S 4 fl Us onto {( G | arg£ = a} and 

{C G /(£/<s) I argC = -cr}, respectively. 



Proposition 5.3 The Riemann-Hilbert problem for P is solved by 

P(z) = E(z)^ a (n 2/3 f(z))e n ^ z)r73 



(5.12) 



1 7 



where 



E(z) = y/ne e 



1 -1 
—i —i 



n 



V*f(z) 



1/4 \ CT 3 



a(z) 



(5.13) 



and ^ a is an explicit matrix valued function built out of the Airy function Ai and its derivative 
Ai' as follows 



* a (0 = { 





( Ai(C) 




V Ai'(C) 




f Ai(C) 




V Ai'(C) 


< 


f Ai(C) 




V Ai'(C) 




f Ai(C) 




V Ai'(C) 



Ai(c 2 C) 
w 2 AiVC) 

Ai(^ 2 C) 
w 2 AiV0 

-w 2 Ai(^C) 
-Ai'K) 

-w 2 Ai(cuC) 
-Ai'K) 



e e 



e e 



<T3 



(73 



1 
-1 1 

1 
1 1 



for < arg ( < a, 
for a < arg £ < 7T, 
/or — 7T < argC < —a, 
/or — a < argC < 0, 



(5.14) 



with u = e 27U//3 . 



Proof. The proof is similar to [10], p. 1523-1525]. 



□ 



Remark 5.4 An analysis of the proof in [TD| shows that the constant C in ( |5.9| ) can be 
taken uniformly for a in a compact subset K a of (vr/3, it), for A in a compact subset Ka of 
(1, oo), and for S uniformly in some interval (So, Si) depending on K a and Ka- 
A similar remark applies to the parametrix P to be constructed below. 

5.3 Parametrix near (3 

A similar construction yields a parametrix P in a neighborhood Us = {z \ \z — (3\ < S} that 
satisfies the following Riemann-Hilbert problem. 

Riemann Hilbert problem for P: 

The problem is to determine, for a given S > sufficiently small, a 2 x 2 matrix valued 
function P : U 5 \ S T -> C 2x2 such that 

(a) P(z) is analytic for z E U$\ S T , and continuous on S T , 

(b) P(z) possesses continuous boundary values for z G S T fl C/j, denoted by P+Oz) and 
P_(V), and 



P+(z) = Pjz) 



1 
-1 



Pjz) = Pjz) 



and 



1 



^ g-2n</>(z) 

1 



p + (^) = p_(^ 

Recall that the 0-function is defined in ( |3.11|) . 



for z g r n c/j, 



for z G (S3 U S 4 ) n Us, 



for z G Si fl U s . 



(5.15) 
(5.16) 

(5.17) 



is 



(c) There exists a constant C > such that for every z G dU& \ £ 



T 



■ 



\\P(z)N-\z)-I\\<-. (5.18) 

n 

There is a one-to-one analytic mapping ( = f(z) such that 

\(-f(z)f"=~4>(z). (5-19) 

Then / maps Us onto a convex neighborhood of £ = 0. 
Proposition 5.5 TTie Riemann-Hilbert problem for P is solved by 

P{z) = E{z)^ a {n 2/3 f{z))e n ^ z)a3 (5.20) 

with 

E(z) = ^e^() l \(n 1 l\-f{z)) 1 / A a{z)Y\ (5.21) 



i —i 



and 

no = ( o A ) ( i -°i ) < 5 - 22 > 

Proof. This follows as in 0, p. 1527]. □ 

Remark 5.6 As in jnj there is a full asymptotic expansion of PN~ l and PN~ l in inverse 
powers of n. This expansion leads to uniform asymptotic expansions for the generalized 
Laguerre polynomials. In the present paper we will compute asymptotics for the polynomials 
to first order, but we will not carry out the further computations to determine the coefficients 
of a complete asymptotic expansion. 

6 Final transformation T \—> S 

Using N, P, and P, we define for every neN, 

S(z) =T(z)N(z)- 1 for z G C \ (S T U Z7 5 U Us), (6.1) 

S(z) = T(z)P(z)- 1 for zeU s \ S T , (6.2) 

S(z) = T(z)P(z)- 1 for z G U 5 \ S T . (6.3) 

Then S is defined and analytic on C \ U dU$ U dUgj . However it follows from the 

construction that S has no jumps on T and on E T PI (Us U Us). Therefore S has an analytic 
continuation to C \ S 5 (also denoted by S), where S 5 is the contour indicated in Figure || 
for A = 1.01. _ 

Formally, the contours are given by = E 3 - \ (Us U for j = 1,2,3,4, and the 

domains flj are given by nf = O x \ (F 5 U f/ 5 ), fif = (H 2 U T U 3 ) \ (Z7 5 U f/ 5 ), and 

n§ = n A \(Usuu 5 ). 

Then S 1 satisfies the following Riemann-Hilbert problem. 



1Q 



-0.5 



0.5 



1.5 



Figure 8: Contour = dll 5 U dU s U [J . and the domains U s , Us, and Qj, j 
the Riemann-Hilbert problem for S, for the value A = 1.01. 



1,2,3, for 



Riemann-Hilbert problem for S: 

The problem is to determine 5* : C \ X 5 — > C 2x2 such that the following hold. 

(a) S^) is analytic for z G C \ S 5 , 

(b) ,S(z) possesses continuous boundary values for z G S 5 , denoted by S+W and S-(z) 
and 



= 5_(z)P(^)iV(2)- 1 for^Gdf^ 



(6.4) 



S + ( z ) = S_{z)P{z)N{z) 



-i 



for z G SiZ? 



(6.5) 



S + (z)=S_(z)iV(z)( \ \N(z) 



-i 



for z G S3 U Xf 



■1 • 



(6.6) 



1 2n0(z) 

S + (z) = S_(z)iV(z) ( q x jiV(z)- 1 for^GSf, 



(6.7) 



S + ( z ) = S-(z)N(z) 



(c) S(z) = 1 + 0(1) for 2^00. 



g -2n0(z) 







x j iV(^)- 1 for G Ef , 



(6i 



The jump matrices for S are close to the identity matrix if n is large. Indeed, by 
and ( |5.18p we have 



\\P(z)N~\z) - I\\ < 



C 



n 



for z G dUs 



9n 



and 

\\P(z)N- l (z) - III < - for z G dU 5 , 

n 

with a constant C that is independent of z. The constant can also be chosen independently 
of the value of A for A in a compact subset of (1, oo), see Remark 5.4. 

The jump matrices in ( |6.6|) - (|6.8|) are exponentially close to the identity matrix. 



7 Asymptotic for polynomials L [ n n) (nz 
7.1 Asymptotics for S 

In Sections 3-6 the values of n and a were assumed to be fixed. In order to study the 
asymptotics, we now let a = a n depend on n and we write A n = —2a. We assume that 
A n > 1 for every n, and 

lim = lim A n = A > 1. (7.1) 

n^oo fl n— »oo 

Because of the n-dependence of A n , all of the notions and results introduced in Sections 3-6 
are n-dependent. For example, we have that the curves T, E, E T , and E 5 are all varying with 
n, and so we denote them by T n , E n , E^, and E^. Likewise, we have that the functions R, g, 
4>, and 0, as well as all matrix- valued functions are n-dependent, and we also use a subscript 
n to denote their dependence on n. If now we use T, E, R, g, etc., without subscript n, then 
this refers to the limiting case. Due to (|7.1|) we have that the curves T n tend to the limiting 



curve T, etc. 

At the end of the previous section, we observed that the jump matrix for S n is l + 0(l/n) 
uniformly on E^ as n — > oo. In addition, the jump matrix converges to the identity matrix 
as z — » oo along the unbounded components of E^ sufficiently fast, so that the jump matrix 
is also close to I in the L 2 -sense. Since the contours Ef are only slightly varying with n, we 
may follow arguments as those given in [0 and ||, to conclude that 

S n (z)=I + 0[-) asn^oo (7.2) 



uniformly for z G C \ E^. Briefly, the Riemann-Hilbert problem for S n is equivalent to a 
system of singular integral equations. This system is of the form 7(1^) — F = g{(), where 
7 is a matrix valued function defined on the contour Ef, and where the singular integral 
operator F(j)(^) has operator norm O(-). The system of singular integral equations can 
therefore be solved by Neumann series in powers of and this leads to a series expansion 
for the solution of the Riemann-Hilbert problem. 

Using more precise information on the jump matrix for S n as indicated in Remark 5.6, 
and with additional assumptions on the limit ( [T.l| ), one is able to obtain a full asymptotic 
expansion for S n (z) in powers of 1/n. This in return would give a full asymptotic expansion 
for the generalized Laguerre polynomials. 



7.2 Strong asymptotics for generalized Laguerre polynomials 

Unraveling the steps Y n \— > U n ^ T n \— > S n and using (|7.2|) , we obtain strong asymptotics 
for Y n in all regions of the complex plane. In particular we are interested in the (1,1) entry 
of Y n , since this is the monic generalized Laguerre polynomial. 
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Theorem 7.1 Suppose for each n e N, we have a parameter a n < —n such that ( |7. 1| ) holds 
where A n = —2a. T/ien we have the following asymptotic results for the generalized Laguerre 



n 

polynomials L^ 1 {nz) as n — > oo. 

(a) (asymptotics away from T) 

Uniformly for z in compact subsets of C\T, we have as n — > oo ; 



-n) 



nz) 



?ng n (z) 



l + R'Jz) 



1/2 



1 + ( - 

n 



(7.3) 



(b) (asymptotics on +-side of T n , away from endpoints) 

Uniformly for z on the +-side of T n away from (3 and f3, we have as n 



nz) 



-n 



ni 



o n 9n(z) 



1/2 



1/2 



3 2n</>„(z) 



+ - 



n 



oo, 



(7.4) 



(c) (asymptotics on — side of T n , away from endpoints) 

Uniformly for z on the —-side ofT n away from f3 and (3, we have as n 



-n 



o n 9n{z) 



/■?■! 



1/2 



X 



l + R'Jz) 



1/2 



e 2n<j) n {z) + Q I _ 

n 



oo, 



(7.5) 



(d) (asymptotics near 0) 

Uniformly for z in a (small) neighborhood of j3, we have as n — > oo ; 



n2 



-n 



exp ^(^logz + z + ^J 



x 



Hi) V4 (^ /3 /n(-)) 1/4 Ai(n 2 / 3 /n(^)) (l + O Q 



• (7-6) 



Proof, (a) Let K be a compact subset of C \ T. Then, for n large enough, say n > n , 
we have i(cC \ T n . Let n > n . While choosing the contours (S 3 ) n and (S 4 ) n in Section 
4, we then take care that (f^)™ and (0 3 )„ do not meet the compact K, see Figure ^|. We 
also choose 5 so that the disks of radius 5 around (3 n and (5 n are disjoint from K. Then 
K C (fif)„ U (fif) n , see Figure g. For z e AT, we then have by Q), Q), and Q 



l n{l n /2)azjj( z \ e ng n {z)a3 e -n{l n /2)(T3 



e n(£ n /2)a3rpr z \ e ng n (z)a3 e -n(e n /2)c 

e n ^ l2)a ' A S n {z)N n (z)e n9n{z)rJz e- n ^ ,2)c 



(7.7) 



99 



Using (|2.8|), (|5.5|), and (|7.2|), and the fact that |(iV n )ii(^)| is bounded away from zero on K, 
we obtain (|7T3| ) from ( |7.7|) . 

(b) For z G (fi 2 ) n \ (([/*)„ U (U s ) n ), we have by (0) and (pj) 

(K„)ii(z) = (tf„)ii(z)e»*»« = [(T„)n(z) + (T n ) 12 (^)e 2 "^] e^M (7.8) 

Since T n = 5 n iV n by (|bM|), = J + 0(l/n), and the entries of N n are uniformly bounded 
away from zero in the region under consideration, we then get 

(y n )nO) = [(AU 11(2) + (AU 12 We 2n<M2) + 0(l/n)] e n ^ z \ (7.9) 

uniformly for z G (0 2 ) n \ ((^)n U (^j)n)- This proves ([f!4[) . 

(c) This is proved similarly as part (b). The only difference is that for z G (fl3)n\((%)nU 
(^<y)n) we use (|4.10p instead of (|4.9|) . This leads to 

(y n ) n (z) = [(T B ) u (z) - (T n ) 12 (^)e 2 "^] e«*»M 

instead of (|7.8| ). The rest of the proof is the same. 



(d) The proof is as in |L0|, p. 1539]. Let U e be an e-neighborhood of j3 with e < 5. Then 
U e C (£/<s) n for all large enough n, say n > n\. Let z E U € and n > Wi. For z G (Oi) n , we 
find by (|OD , §3), and (gj) that 

{Y n ) n (z) = (U n ) xl {z)e^ = (T n ) u (z)e n ^ 

= [(S n ) n (.z)(P n ) n (z) + (S n ) 12 (z)(P n ) 21 (z)] e n9 '^. (7.10) 
By Proposition 5.3 we have 

P n (z) = E n (z)^(C)e n ^ z ^ 

where £ = n 2 ^ 3 f n (z). For z G (£)i) n , we have that £ belongs to the sector < arg£ < er, so 
that the first formula for \I/ CT (C) in Q5.14Q applies. It follows that 

(P n )n(z) \ _ En{z) f Ai(C) \ e -f e nM*), (7 .ii) 



(Pn)u(z) J ~ V Ai'(C) 

Using the definition ( |5.13| ) of E n (z), we obtain from ( [7 . 1 1| ) 

V {Pn)x2{z) ) Vn { -i -1) {a n (z)J \ Ai'(C) 



Then 



{S n ) u {z){P n ) n {z) + {S n ) l2 {z){P n ) 21 {z) 



4^Ai(o 

a„(z) \W 



^(1 + 0(1) -1+0(1) ) ( |e"*»W. (7.12) 

$?Ai'(C) 

Combining ( |7A0| ) and ( [7l2|) with ( = n 2 l 3 f n (z), we obtain an expression for (Y n ) u (z), which 
leads to Q in view of Q and (pT3|) . 

Next, for z in the other regions (%) n > J — 2, 3, 4, similar calculations lead to the same 
expression for (Yu) n (z). Hence ( \l.6\) holds uniformly for z G U e . 

This completes the proof of Theorem 7.1. □ 



Remark 7.2 Since a n is real, we have L^ n \nz) = Ll? n \nz), and so (|7.6| ) also describes 
the asymptotic behavior near /3. 
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7.3 Asymptotics for zeros 



From Theorem 7.1 we deduce the following results concerning the zeros of the generalized 
Laguerre polynomials L^ n \nz). 

Corollary 7.3 Suppose we are in the same situation as in Theorem 7.1. 

(a) (All zeros tend to T) For every neighborhood Q ofT, there is no such that for every 
n > n , all zeros of L^ n \nz) are in Q. 

(b) (Zeros are on the — side) For every 5 > 0, there is no such that for every n > no, 
there are no zeros of Ln (nz) in the region (f^n \ ((Us) U (Us)). 

Proof, (a) This is immediate from the asymptotic formula ( |7.3| ) since 1 + R' n (z) ^ 0. 

(b) Suppose z is a zero of L^ n \nz) lying in (fi 2 ) n \ ((Us) U (Us)), that is, on the +-side 
of T n . Then we have by (|7.4j) 

( i— V2 e 2nMz) = l + o(-). (7.13) 
\l+R' n (z)J \n) y } 

We know that Re(j) n (z) < 0. In order to obtain a contradiction it is thus enough to enough 
to prove that 

\1 - R' n (z)\ < \1 + R' n (z)\. (7.14) 

Equality holds in (|7.14|) if and only if R' n (z) is purely imaginary, so that (R' n ) 2 (z) is real 
and negative. From the definition ( |3.2| ) of R n , we get 

y^x (z-^f (z-ReP n f 



(KY(z) 



(z - (3 n )(z - (3 n ) (z - Re(3 n y + (lm/3 n y 



This is real and negative if and only if z belongs to the vertical segment connecting (3 n and 
(3 n . Consequently, this is the set where \1 — R' n (z)\ = |1 + R' n (z)\. The vertical segment 
and T form a closed contour, and it may be checked that ( |7. 14|) holds for z outside of this 
contour, which includes the region (f2 2 ) n \ ((Us) U (Us)). 

So we have a contradiction, and it follows that there are no zeros in (f2 2 )n \ ((Us) U (Us)) 
for n large enough. □ 



Remark 7.4 From the uniform asymptotics ( |7.6| ) in a neighborhood of (3 it is possible to 
derive asymptotics for the extreme zeros of L^n n \nz). For example, it follows that for fixed 
k G N and n > k, there is a zero z^^ n of Ln (nz) such that 

^ = A --fi^ + ©- (7 ' 15) 

where —tk is the fcth largest (negative) zero of the Airy function Ai(x). 

OA 
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